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Abstract 

In this paper, we establish a set of criteria which are applied to discuss various formulations 
under which Lagrangian stochastic models can be found. These models are used for the simulation 
of fluid particles in single-phase turbulence as well as for the fluid seen by discrete particles in 
dispersed turbulent two-phase flows. The purpose of the present work is to provide guidelines, 
useful for experts and non-experts alike, which are shown to be helpful to clarify issues related to 
the form of Lagrangian stochastic models. A central issue is to put forward reliable requirements 
which must be met by Lagrangian stochastic models and a new element brought by the present 
analysis is to address the single- and two-phase flow situations from a unihed point of view. For that 
purpose, we consider first the single-phase flow case and check whether models are fully consistent 
with the structure of the Reynolds-stress models. In the two-phase flow situation, coming up 
with clear-cut criteria is more difficult and the present choice is to require that the single-phase 
situation be well-retrieved in the fluid-limit case, elementary predictive abilities be respected and 
that some simple statistical features of homogeneous fluid turbulence be correctly reproduced. This 
analysis does not address the question of the relative predictive capacities of different models but 
concentrates on their formulation since advantages and disadvantages of different formulations are 
not always clear. Indeed, hidden in the changes from one structure to another are some possible 
pitfalls which can lead to flaws in the construction of practical models and to physically-unsound 
numerical calculations. A first interest of the present approach is illustrated by considering some 
models proposed in the literature and by showing that these criteria help to assess whether these 
Lagrangian stochastic models can be regarded as acceptable descriptions. A second interest is to 
indicate how future developments can be safely built, which is also relevant for stochastic subgrid 
models for particle-laden flows in the context of Large Eddy Simulations. 
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I. INTRODUCTION 


Over the last decades, Lagrangian stochastic models have become increasingly used for 
both single-phase reactive flows and dispersed two-phase turbulent flows (with one phase 


de^ . These ap- 
which indi¬ 


being present as discrete elements such as solid particles, droplets or bub 
proaches are referred to as PDF (Probability Density Function) methods 
cates that they are simulating the pdf of the relevant variables which have been retained 
for the statistical description of either single- or two-phase flows. PDF methods have strong 
advantages as they treat important phenomena without approximationdn single-phase reac¬ 
tive flows, this corresponds to convective and reactive source terms l|, l9|] while, for dispersed 
two-phase flows (even for inert particles), this corresponds to tran^ort and polydispersity 


effects (related to the existence of a range of particle diameters) [l^ . This interplay of 


modeling issues explains the common interest of PDF descriptions in both situations. 

In the single-phase flow situation, the governing equations are the transport equations 
for the fluid velocity held U(t,x) and for a set of scalars which gathers the relevant species 
mass fractions (/>(t, x) = ( 0 / 3 )/ 3 =i,...,Ar^ to which an equation for the huid enthalphy is added, 
along with an equation of state, for compressible hows. For constant-property hows, these 
equations are 


dUk 

dxk 


= 0 




dUi dUi 



dt 



IdP ^ d‘^Ui 

pdxi dxkdxk 


dxkdxk 


+ S0 


(la) 

(lb) 

(l c) 


where u is the huid dynamical viscosity, F the scalar dihusivity and P(t, x) the huid pressure. 
In Eq. flTcl) . the last term on the rhs (right-hand side) is the reactive source term = 
Sj3{4>{t,x.)) which, along with convection, appears in a closed form in a one-point PDF 
approach. 

In the disperse two-phase how case, the basic physical situation is made up by a continuous 
huid phase (a gas or a liquid) in which a set of discrete ‘particles’ (solid particles, droplets, 
bubbles, etc.), having a range of diameters, are embedded. The huid phase is described by 
the continuity and Navier-Stokes equations, Eq. (ITaD -Eq. flTbl) . to which source terms can 
be added to account for momentum exchange between the huid and the discrete particles 
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when two-way coupling (whereby particles influence the fluid phase) is deemed important. 
For the discrete particles, we limit ourselves to the case of point-like particles or droplets. 
This approximation is not severe for most industrial applications but usually leaves out 
bubbles [n|. For particle diameters of the same order of magnitude as the Kolmogorov 
length scale, the particle momentum equation involves the well-known pressure-gradient. 


drag, added-mass and Basset forces [l2|, ll3| . In the case of particles heavier than the fluid 
(droplets in a gas, solid particles in a gas or liquid), the particle momentum equation can 
be simplified to the following form to describe the evolution of particle location Xp(t) and 
velocity Up(t): 

dyir, 


dt 

dt 


= Up, 

- = -(U, - Up) + g. 


Tp 


(2a) 

(2b) 


where, apart from gravity, only the drag force has been retained (other forces can be added 
but the drag force is sufficient for the present discussion). In Eq. (l2b|) . g is the gravity 
acceleration and Tp the particle relaxation time defined as 


Tp = 


Pp 


4:dri 


( 3 ) 


Pf SUi^lUrl ’ 

where the local instantaneous relative velocity between the fluid and the particle velocity is 
Ur = Us — Up. The drag coefficient Ud is usually expressed as a non-linear function of the 
particle-based Reynolds number, Rcp = dp\\Jr\/i^ (where dp is the particle diameter) based 
on empirical formulas Q (apart from the Stokes regime). In Eq. fl2bj) and in the expression 
of the particle relaxation time scale, the important variable is Us(t) = U(t,Xp(f)) which is 
the ‘fluid velocity seen’, i.e. the fluid velocity sampled along the particle trajectory Xp(t) as 
it moves across a turbulent flow. In the limit of the assumptions made above, it is possible to 
solve the governing equations in the spirit of DNS (Direct Numerical Simulation) for disperse 


flows with accurate Lagrangian tracking methods ISMlSl]. as the continuation of what is done 


for the DNS of single-phase flows. This approach is possible in simple geometries and at 
moderate Reynolds numbers. When such DNS are not available and for practical purposes 
where only one-point statistics are known, a stochastic model for the velocity of the fluid 
seen is needed. In the general context of stochastic models and PDF descriptions, we are 
thus concerned with a one-particle PDF model for Us(t) from which one-point statistics can 
be derived jd]. Note that when the velocity of the fluid seen is included in the state-vector 
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along with the particle velocity and diameter, the drag force appears in a closed form in the 
PDF description jdj (it can be seen as the counterpart of the reactive source term in the 
single-phase flow case). 

In the present analysis, we leave out the specihc issues related to reactive aspects (gas 
chemical reactions in single-phase flows; burning particles, evaporating droplets, etc. in two- 
phase flows) and concentrate on dynamical aspects in the PDF approach to both situations. 
Thus, at the core of the PDF method lies the specihc stochastic model which is used to sim¬ 
ulate huid particle velocities in single-phase hows and the huid velocity ‘seen’ by particles in 
two-phase hows. In the latter situation, this notion was introduced as a Lagrangian property 
attached to each discrete particle 1^, then included in the particle state-vector 10| which 
led to the standard PDF description for two-phase hows Furthermore, the PDF ap¬ 


proach was given a complete framework, hrst in single-phase hows [l| and later in two-phase 
hows jd], ensuring a continuous link between the choice of the particle state-vector, the for¬ 


mulation of stochastic models as proper stochastic diherential equation^ 
PDF equation in sample space and the resulting mean-held equations [1, 


;he corresponding 


20 


2l|. 


At this stage, a hrst modeling question arises: how should these stochastic models be 
formulated? In single-phase hows, most of tl^ stochastic models have been developed in 
terms of instantaneous huid particle velocities 


20 


22| . In two-phase hows, some proposals 


followed the same road and were made in terms of the instantaneous huid velocity seen to 
build so-called Langevin models jd, 23|. However, the modeling situation is unclear since 


existing models can be expressed with diherent formulations, for instance in terms of either 
instantaneous, huctuating or normalised huctuating huid-particle velocities. When particle 
inertia becomes negligible, the stochastic model for the velocity of the huid seen becomes 
a stochastic model for huid particles and, consequently, this issue overlaps with similar 
concerns in single-phase hows. This means that both situations are impacted by the issue 
of the formulation of stochastic models. A second modeling question is: what are the basic 
properties that such models must respect? 

These questions are particularly relevant when the stochastic models used to simulate 
dynamical variables are complemented with additional models to address complex-physics 
problems. For instance, in single-phase hows, one could be interested in applying a velocity- 
composition PDF method but having only specihcally developed the modeling parts con¬ 
cerned with scalars and reactive terms. In two-phase hows, Lagrangian stochastic models 
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can be used to analyze additional effects such as thermophoresis, electrophoresis or chemical 


forces 


24 


(due to interface chemistry in liquid medium). Other typical examples include 


droplets or coal/fuel particles where models are added to simulate complex combustion or 


evaporation processes 


flis, 


m 


Yet, these practical developments can be ruined by a 


poor formulation of the model retained for the velocity of fluid particles or the velocity of 
the fluid seen. It appears therefore important to assess whether stochastic models used for 
dynamical variables have a sound basis. By this assessment step, we do not mean here the 
(necessary) task of outlining the predictive abilities of modeling proposals but the ‘upstream 
step’ that consists in assessing whether the structure of these models respect key properties, 
regardless of the details of specihc closures. Going directly to a comparison of computa¬ 


tional outcomes between different model formulations 


29| is interesting but may confuse 


the status of different classes of Lagrangian stochastic models whereas a theoretical analysis 
can already reveal flaws or bring out differences that make these models difficult to compare 
even-handedly. 

Given the subtleties of stochastic calculus 
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3l| and the particle-based nature of La¬ 


grangian PDF approaches, consistency issues have accompanied the construction of stochas¬ 
tic models. For example, the issue of so-called spurious drifts and its relation to fluid mass 
conservation was addressed a few years ago for models developed for single-phase turbulent 
Relations with Reynolds-stress modeling were also put forward 20| with a 


flows 


32 


view towards the interest of stochastic mode 
been made in the two-phase flow situation 


s for realizable closures. Similar efforts have 
, though more sparingly. However, there has 


been no previous attempt at gathering knowledge and addressing the validity and the struc¬ 
ture of stochastic models by resorting to a systematic list of requirements, especially when 
the single- and two-phase situations need to be jointly considered. 

With respect to this context, the hrst purpose of this article is to propose a clear set 
of requirements, for single-phase as well as for two-phase flows, which must be met by 
Lagrangian stochastic models in order to be regarded as acceptable descriptions. The second 
purpose is to discuss the relations between different formulations of a stochastic model for 
the velocity of fluid particles and to reveal the interests as well as the limitations of some 
of these formulations. The third purpose is to establish guidelines for future developments, 
valid for stochastic models developed in classical Reynolds-averaging approaches but also 
of interest for models considered for subgrid-effects in particle-laden flows where the fluid 
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phase is simulated with a LES (Large Eddy Simulation) approach. In that sense, the present 
considerations represent an effort to address issues related to single-phase and two-phase 
ffows from a unified perspective. 

The paper is organized as follows. The PDF theoretical framework is first recalled in 
section m where the simplifying assumptions which define the precise context of this study 
are stated at the end of section IIIBl Then, the criteria selected in the present analysis are 
detailed in section HU first for the single-phase flow case in section UlI Al and, second, for the 
two-phase flow case in section UlIBI The analysis is first developed for single-phase turbulent 
ffows: different formulations for fluid particle velocities are addressed in section IIVI and 
analyzed in detail in section IIV Al while classical scalar modeling (which plays an important 
role in one criterion for two-phase ffows) is recalled in section IIVBI Drawing on the analysis 
carried out for the single-phase flow case, an analysis of different modeling proposals for two- 
phase ffows is developed in section |Vl In particular, new relations for models expressed in 
terms of fluctuating components, as well as discussions on two-way coupling, are developed 
in section IV A 31 Present findings are summarized in Table [I] for single-phase models and in 
Table m for two-phase models. Finally, guidelines for future developments are proposed in 
the Conclusion. 


II. THE PDF THEORETICAL FRAMEWORK 


This framework was first established for single-phase turbulent ffows [l-^ 


20 and was 
vs 3,21,36|. 


later used as a foundation for the extension to dispersed two-phase turbulent ffows 
However, for the sake of a simpler and more compact presentation, we introduce the key 
aspects of the theoretical framework directly from the standpoint of the two-phase flow 
situation since the fluid-particle case can be retrieved as an asymptotic limit. 


A. Probabilistic descriptions and stochastic equations 

The PDF machinery for fluid mechanics starts by the choice of the PDF description 
(in terms of either one-particle pdf, or two-particle pdf, etc.). In the present context, we 
consider only one-particle PDF approaches and, thus, the starting point is the selection of 
the mechanical description retained for each particle or, in other words, the choice of the 


7 























relevant particle state-vector which gathers the variables of interest attached to each particle. 
Following the presentation of standard Lagrangian models for dispersed two-phase flows in 
the Introdnction, the particle state-vector is made np by the particle location and velocity 
as well as the velocity of the flnid seen by the particle, Z = (xp, Up,Us), with evolntion 
eqnations written as 


dxp = Up dt, 

dUp = Dp(t, Z) dt, 

d\Js = D,(t, Z, ^[(Z)], ($)) dt + B,(t, Z, ^[(Z)], ($)) dW. 


(4a) 

(4b) 

(4c) 


In these eqnations, Dp typically represents the drag and gravity forces. Dp = (U<j — Up) /Tp-f 
g (other forces can also be considered), while the vector D^ and the matrix are the drift 
and diffnsion coefficients of a stochastic diffusion process which is a typical model for the 
velocity of the fluid seen. In the drift and diffusion coefficients in Eq. fl4c]) . a general notation 
has been used to indicate that these coefficients can depend on the value of the state-vector 
Z but also on functionals of the mean fields which are calculated from the simulation of that 
state-vector, written as J^[(Z)], as well as on some external helds represented by (<F). A 
typical example of J^[(Z)] is the particle mean-velocity held while the huid mean-pressure is 
another example of what (<h) can stand for. For detailed discussions of the modeling issues 
from a physical point of view, we refer to the relevant literature { 2 , 4|. 

The single-phase how framework is retrieved by considering the particle-tracer limit when 
particle inertia goes to zero (up —0). In that case, the particle velocity tends towards the 
huid velocity and, for example, the model system of equations, Eqs. (0]), becomes 


dx = U dt, 

dV = D(t, Z, -F[(Z)], ($)) dt + B{t, Z, -F[(Z)], ($)) dW. 


(5a) 

(5b) 


where the same notation has been retained to indicate a possible dependence on mean-helds 
calculated from the solution (such as the huid mean velocity held) or on external mean 
helds. A typical example of such an external mean held is the huid mean dissipation held, 
(e). It must be noted that the notion of external helds is, of course, directly dependent 
upon the choice of the variables entering the state-vector. For instance, for Z = (x, U) 
which corresponds to a velocity-PDF description 22|, the mean dissipation is an external 


held. Yet, if the state-vector is extended to include the instantaneous dissipation, whereby 






Z becomes Z = (x, U, e) and corresponds to a velocity-dissipation PDF description [91,1371, 


0 , 


then the dependence of the drift and diffusion coefficients on the mean dissipation held would 
appear through the functional form J^[(Z)]. 

For both the single- and two-phase situations, using a general equation written 


dZ = 'Dzdt + Bz dW 


( 6 ) 


as a reference model, the key points of the PDF theoretical framework can be unraveled. 

The governing equation, Eq. ([6]), is to be understood as the evolution equation for a 
large number of ‘stochastic particles’ (whose behavior mimic the evolution of real particles 
in a statistical sense) and, in a weak sense [l, 4, 38, 39|, a stochastic particle-tracking 


model is equivalent to a PDF approach. However, it is worth emphasizing that the present 
case is an extension of the well-established connection between classical Langevin equations, 
where the drift and diffusion coefficients depend only the chosen state-vector Z, and the 
corresponding Fokker-Planck equation in sample space. In that sense, there is an essential 
difference between the Langevin equations used, for example, in PDF methods and in dis¬ 
persion studies 0. The difference is that, in PDF methods, the coefficients involve statistics 
obtained from the particles as indicated by the general notation in Eqs. (jl]) and ([5]). In the 
mathematical literature, these processes are referred to as ‘McKean diffusion processes’ 
and in detailed physical presentations ^ they are defined as ‘processes with mean-field 
interactions’. Compared to classical Monte Carlo methods, numerical implementations lead 
to handling so-called ‘weakly interacting processes’ but, basically, the classical connection 
between these generalized Langevin and Fokker-Planck equations remains valid {m, 36|. In 
other words, a particle-tracking method amounts to a dynamical Monte Carlo simulation of 
the corresponding pdf pL(f;z), where the index L indicates that we are dealing with a La- 
grangian pdf and where z stands for the state-vector variables in the corresponding sample 
space. 


B. Prom Lagrangian stochastic models to mean-field equations 

In order to go from the (Lagrangian) stochastic equations written for each notional parti¬ 
cle to the resulting (Eulerian) mean-field equations, a straightforward approach is obtained 
by introducing MDFs (Mass Density Function) 0 and by dehning the Eulerian MDF from 
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the Lagrangian one. This is translated in the following dehnition where the Eulerian MDF is 
basically identihed with the Lagrangian one at the same loca^on (the identity is for MDFs 
and not necessarily in terms of pdfs) as detailed elsewhere 
mass of the physical particles in the domain: 


^ and where Mp is the total 


Fp{t;yp,Zc) = MpPL{t;yp,Zc), ( 7 ) 

Fp{t, x; zj = Fjf (t; = x, zj = j yp, z^)6{yp - x) dy. (8) 

Since it is useful to distinguish between the particle location and other variables, we have 
introduced the notation Z = (xp, Zc) for the particle state-vector where the particle location 
is always present and where Z,, stands for the complementary part of the chosen state-vector. 
For example, for the Lagrangian description mentioned above for dispersed two-phase flows, 
we have Zc = (Up, Us), while for the single-phase flows we would have Zc = (U). In 
sample-space, the corresponding variables of the state-vector are noted z = (yp,Zc). 

For a particle variable written as iLp(t; Zc), its average (iLp)m (which is a held variable), 
is dehned as 

ap{t, x) pp{Hp)m{t, ^) = j '^c)Fp{t, x; Zc) dz^. (9) 

where Pp is the particle density and ap{t,x.) is the mean particle volumetric fraction. The 
huctuating component is then expressed as hp = Hp — {Hp)m- The dehnition of these 
quantities for the single-phase case are given below, after Eq. (ITT)) . It is worth emphasizing 
that Op is a rigorously-dehned probabilistic quantity that represents the average presence 
of one phase at a given location and should not be confused, in the present framework, 
with volumetric averages. In a discrete sense, when we handle N stochastic particles, the 
dehnitions of Lagrangian and Eulerian MDFs are directly carried out to yield 

N 

Fp,Nit ; yp, Zc) = 5^mW5(yp - x^) 5(zc - Z^) (10) 

i=l 

Fp,Nit, X; Zc) = Fp^^(t; yp = X, Zc) (11) 


where rrip'^ is the mass of the particle labeled {i). This shows that in a small volume around 
location x where averages are estimated as the ensemble averages over the particles 
present in that volume, we get the equivalent of Favre, or mass-weighted, averages 


{Hp)^{t,x.) ~ (Hp) 






m,N 


Z^i=l 


(0 

my 


( 12 ) 
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which explains the notation {Hp)m used to indicate mass-weighted averaging for particle 
variables. 

One of the key points is that the Eulerian MDF satisfies the same evolution equation as 
the Lagrangian MDF. For the general stochastic diffusion model considered in Eq. dH]), this 
equation is a Fokker-Planck equation 


m dzk ] + 2dztdzi ] ■ 


(13) 


For the standard state-vector for dispersed two-phase flows, Z = (xp,Up,Us), and for the 
evolution equations considered in Eqs. dl]), the Fokker-Planck equation has the following 
form where convection appears in closed form: 

(14) 


dF^ dF^ 


dt 


p,k' 


dyk 


dVp^k 

d 


dK 


s,k 


(D.m Ff) + 


8^ 


2dV,_tdV.i 




where the variables in sample-space are denoted (y, Vp, V^). It is instructive to consider the 
first particle-velocity moments obtained from the PDF methodology 


doipPp 0 (^OipPp{Up^k')m) p, 

d{Up,)r 

dt 


dxk 
4 “ {Up,k')i 


d{Up,^)r 

dxi 


1 d {otpPpiUp^iUp^l^ 
dxi 


+ 


U^.i - D, 


p,l 


+ 9i 


(15) 

(16) 


j]^ OLpPp ^'-^k \ ‘p 

where {up^iUp^k)m is the particle kinetic tensor. In the rhs of the average particle momentum 
equation, Eq. flT6|) . the drag and gravity force terms have been directly introduced instead 
of the more general notation Dp. This helps to reveal that in two-phase flows the closure of 
the drag force term is akin to the closure of the chemical source term in single-phase reactive 
flows. Indeed, even for inert particles, Tp is a complex function of the particle diameter and 
of the particle and fluid seen velocities (cf. Eq. d3])). Thus, for a polydisperse set of particles, 
the last term on the rhs of Eq. ffT6|) cannot be closed directly whereas it is handled without 
approximation with the PDF approach. It is also interesting to consider the second-order 
equations for the particle kinetic tensor components which correspond to the Reynolds-stress 
equations in single-phase flows. The exact forms of these equations are detailed elsewhere 


21 | and have an intricate structure which implies further closure issues, such as {Dp^iUpj)m, 


apart from classical third-order correlation terms. Even when the set of polydisperse particles 
is separated into classes of monodisperse ones within which Tp can be regarded as a (nearly) 
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constant parameter, it is readily seen that, for each class, the eqnations for {up^iUp^k)m 
involve the (non-symmetrical) flnid-particle velocities correlation tensor {us,iUpj)m whose 
eqnations, in tnrn, rely on the tensor {us^iUsj)m derived from the velocity of the flnid seen. 
The important points are that these complete second-order eqnations (which represents 
typically 25 conpled partial differential eqnations (PDE) for each particle class) are actnally 
derived from the stochastic description and involve both the drift and diffnsion coefficients 
of the model nsed for the velocity of the flnid seen in Eq. flT^ . This indicates that safe 
gnidelines are needed both to assess existing models and to snggest proper formnlations. 

In the single-phase flow sitnation, all the relations presented above remain valid. In 
that case, Pp is eqnal to p the flnid density, Mp is replaced by M/ the total flnid mass 
in the domain, becomes the mass associated to each notional flnid particles and 
Q!p(t, x) = 1 at any point x. Some comments may be added for the last two points. In 
dispersed two-phase flows, the particles handled in the nnmerical simnlations represents 
stochastic particles bnt since they mimic the ‘real’ particles carried by the flnid tnrbnlent 
flow, their diameter and their mass can correspond to the real-particle properties (provided 
no statistical weights are nsed whereby stochastic particles represent so-called parcels or 
gronps of particles). In that sense, the notion of a stochastic particle may be less snrprising 
at hrst sight in the two-phase case than in the single-phase one, althongh the notions are 
basically the same. In the single-phase case, the mass associated to each stochastic particle 
appears more as a parameter representing a fraction of the total mass, for example eqnal 
mass can be assigned to each particle as often done for incompressible flows. Fnrthermore, 
when the flow occnpies the whole domain, the natnral limit is to have ap{t, x) = 1. Actnally, 
it will be seen in the next section that this corresponds to the absence of spnrions drifts. 
With these dehnitions, the above formalism reverts to the developed for single-phase reactive 


flows 


M- 


In the following, a nnmber of simplifying assnmptions are made. First of all, we consider 
incompressible flnid flows, where the flnid density p is constant. In correspondence, particle 
densities are also taken as constant and, fnrthermore, we assnme that all particles have the 
same constant density (the case of particle-dependent densities is natnrally included since the 
formalism already accounts for particle-dependent masses). However, it must be noted that, 
in two-phase flows, variable particle volumetric fractions ap{t,x.) always induce variable- 
density effects, which justihes the introduction of MDFs as done in this section. In the 
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present analysis, we limit ourselves to high Reynolds-number flows where viscous effects can 
be neglected (apart from the non-zero dissipation rate of the fluid turbulent kinetic energy) 
and, similarly. Brownian effects are not considered for particle motion. In the context of 
this first analysis, this is done to limit the complexity of the stochastic models and to 
concentrate on the important convective and dissipative issues. Furthermore, few models 
have been devised for low-Reynolds particle-laden flows. Note, however, that formulations 


in the case of low Re-number flows were proposed for near-wall fluid flow simulations 


41 


43| 


and that extensions of the present analysis can be made. For the two-phase flow situations, 
two-way coupling is not explicitly accounted for, although related aspects are discussed 
when the form of the equations for the fluctuating velocity components are considered (see 
section IV A 3p while the effects of particle-particle collisions are left out. 


III. CRITERIA FOR THE ANALYSIS OF STOCHASTIC MODELS 

We consider Lagrangian stochastic models for fluid particles, such as those characterized 
by Eqs. ([5]). The criteria put forward here arise from the following question: in what sense 
can we regard these models as acceptable probabilistic descriptions of turbulent flows? 


A. Choice of criteria for single-phase flow models 


fl 


This issue was addressed by Pope [l|, l33|] who showed that, for such models to be true 


PDF descriptions (that is for the corresponding pdf to be normalized to one), it is necessary 
and sufficient that the mean-pressure gradient be properly introduced and that the mean 
continuity condition be respected. This is tantamount to saying that Lagrangian stochastic 
models are free from spurious drifts when the mean pressure gradient is properly taken into 
account { 4 , 33, 35|. The correspondence between Lagrangian stochastic models and Reynolds- 


stress closures was then brought forth 


2d, 


22| . Drawing on these relations, we propose to 


regard a Lagrangian stochastic formulation as an acceptable model if three criteria are met: 


(F-1) the model should be capable of being written in correct Cartesian tensor notation 
and should satisfy the relevant invariance principles (such as Galilean and extended 
Galilean invariance); 


13 


















(F-2) the stochastic model should be such that the mean-continuity equation is respected 
and, thus, the model must be free of spurious drifts; 

(F-3) all convective terms in the mean Navier-Stokes and in the second-order equations 
must be exactly reproduced (as well as dissipative terms). 


The hrst criterion amounts to asking that the basic transformation and invariance properties 
of the Navier-Stokes equations be reflected in the model j^. The last two criteria mean that 
the hrst two moments of the velocity helds extracted from the simulated pdf are such that, 
using Reynolds decomposition written as Ui = {Ui) + Ui, we have the following structure 


d{Uk) 

dxu 

dm 

dt 

d{uiUj) 

dt 


= 0 




dxi 


dxi 


l d{P) 
P dxi 




dxk 


dxk 


dxi 


dxi 


(17a) 

(17b) 

(17c) 


The hrst equation is the mass-continuity equation when the huid density is constant. The 
second equation corresponds to the mean Navier-Stokes equation (i.e., the Reynolds equa¬ 
tion) and it is seen that the mean viscous term, uAUi, has been neglected, for the sake 
of simplicity in the present analysis, by assuming sufficiently high Reynolds-number hows. 
This indicates that we are not considering near-wall regions. In Eq. fll7cp . the rhs (right- 
hand side) Tjj is usually decomposed as Tjj = ^ij — eij where and Cij stand for models 
expressing the correlation between the pressure-gradient and huctuating velocities and the 
dissipation sink term, respectively. These models are subject to classical requirements, for 
instance that = d{pui)/dxi or, when the divergence of the pressure-velocity correlation 
is neglected, to <Fjj = 0 and usual model simplihcations, such as isotropic dissipative terms, 
which means that we can write eij = —2/3{e)6ij with (e) the mean turbulent kinetic energy 
dissipation. In the following, we only require that half the trace of the rhs has the form 
1/2 Tjj = V./ — e with e > 0 and where I stands for an expression that represents the 
pressure-fluctuation correlation. As a practical matter, pressure transport is often ignored, 
and so, for simplicity, we take the requirement to be 1/2 Tjj = —e. 

It is worth noting the difference with the previous works that established the correspon¬ 
dence between Lagrangian stochastic models and Reynolds-stress equations [20|, l22|. In 


these works, the aim was mostly to obtain detailed forms of Reynolds-stress models from 
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the Lagrangian stochastic approach and, in that sense, the emphasis was pnt on developing 
detailed closnre proposals for and e^-. In the present work, the emphasis is pnt on a 
different point. Indeed, by the criterion (F-3) listed above, it is meant that all the terms ap¬ 
pearing on the Ihs (left-hand side) of Eqs. flT7)) shonld be exactly reprodnced. This is a rather 
natnral reqnirement to set forth: indeed, all the terms gathered on the Ihs of Eqs. fll7p are 
convective terms which arise from the transport term in the Navier-Stokes eqnations. Since 
Lagrangian approaches are precisely attractive as they treat convective transport withont 
approximation, it is therefore natnral to expect that the strnctnre of the transport eqnations 
for the hrst two moments of the velocity held will be exactly reprodnced. In other words, 
failnre to do so by a Lagrangian model means that transport is poorly described and this 
is regarded as an nnacceptable shortcoming. The criteria proposed here are basic physical 
reqnirements bnt represent a step forward with regard to former proposals 341. 


B. Choice of criteria for two-phase flow models 


We now consider Lagrangian stochastic models for the velocity of the flnid seen by discrete 
particles, as expressed by Eqs. (jl]), and the issne is to come np with relevant criteria to assess 
whether these modeling proposals can be retained as acceptable descriptions. However, 
selecting snch a list of criteria for two-phase flows is not straightforward. Indeed, the analysis 
of single-phase stochastic models can rely on the sonnd basis provided by given Reynolds- 
stress models. This is not so in the two-phase sitnation as it is actnally the stochastic model 
which is nsed to derive corresponding mean-held or so-called continnnm descriptions of two- 


phase flows 


a 
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44l |. This lack of a sonnd reference continnnm theory pnts a stronger 


emphasis on the need to have reliable criteria to assess stochastic models for disperse two- 
phase flows. In spite of this limitation, it is still possible to nse the single-phase flow case 
as a sonnd basis. More specifically, to analyse stochastic models nsed for the simnlation of 
the flnid velocity seen Us, we propose the following list of criteria: 


(P-1) the model shonld be capable of being written in correct Cartesian tensor notation 
and shonld satisfy the relevant invariance principles (snch as Galilean and extended 
Galilean invariance); 

(P-2) the stochastic model mnst be complete in the sense that the expressions for both the 


15 







drift vector and the diffusion matrix in Eq. must be explicitly given. Furthermore, 
the stochastic model must have predictive capacities in the sense that it should be 
applicable to general non-homogeneous situations where fluid or particle statistics are 
not necessarily known in advance; 

(P-3) the stochastic model used for the fluid velocity seen should revert to an acceptable 
stochastic model for single-phase flow (thus respecting the criteria (F-1) to (F-3)) when 
particle inertia goes to zero, that is in the particle-tracer limit; 

(P-4) the complete stochastic model should be such that the predicted mean kinetic energy 
of the fluid seen respects the same statistical evolution laws as the mean fluid kinetic 
energy in dilute flows, both in stationary and decaying homogeneous turbulence; 


(P-5) the resulting model for the so-called drift velocity, say (which represents the dif¬ 
ference between the mean velocity of the fluid seen and the mean fluid velocity at 
the same particle location, or Vd,i = {Us,i)m — {Ui)), should be consistent with known 
scalar dynamics and related modeling; 


(P-6) the model for the velocity of the fluid seen should be consistent with the Equilibrium 


Eulerian approach formulated in the limit of small particle Stokes numbers 18|, |45 |. 


More precisely, the model should be such that the mean conditional increments of the 
drift vector yield the correct behavior for discrete particle velocities for small Tp. 


The criterion (P-1) is the same criteria as for the fluid case (see (F-1)) carried to the 
two-phase flow situation. Indeed, since the velocity of the fluid seen is extracted from the 
velocity held of the huid carrier phase, the same invariance principles hold. At hrst sight, the 
criterion (P-2) seems obvious as it simply express what a ‘model’ stands for. Unfortunately, 
it will be indicated in section |V] that there can be some confusion as to which proposals 
can be truly referred to as a stochastic model or, in other words, as a true PDF description 
of two-phase hows. Thus, this criterion will be shown to be useful to clarify the status of 
some proposals. To be self-consistent, it is useful to clarify the dehnition retained here for a 
predictive model: a model is said to be predictive if it is able to provide information on the 
future state of a system (which means on the statistics of the variables retained to describe 
the two-phase how), given all necessary initial and boundary conditions. 
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In the list above, the criterion (P-3) is an important one: it translates what is meant when 
it is said that satisfactory models for two-phase flows should be extensions of satisfactory ones 
for single-phase flows. From a physical point of view, this is a sound requirement to make: 
indeed, small-inertia particles tend to follow closely the surrounding fluid and it appears 
consistent to expect that statistics derived from the fluid and from the discrete particle set 
become identical and both acceptable descriptions. Furthermore, from a numerical point 


of view when so-called Eulerian/Lagrangian hy 
ensure that consistent predictions are obtained 


Drid simulations are used, this will help to 


46|. 


The last criteria involve specihc aspects of two-phase flows. The criterion (P-4) states 
for a statistically homogeneous decaying turbulence, we expect that, for particles which are 
homogeneously distributed (still in a statistical sense), the fluid seen is such that 

1 d{us ■ 


dt 


= -w ■ 


( 18 ) 


This expresses that the mean fluid kinetic energy seen by particles is decaying at the same 
rate as for the fluid. In the context of the present study, this criterion is meant as a functional 
requirement and it will be seen in section |V] that it is particularly useful to lead to correct 
closure expressions for the diffusion term in Eqs. (0]). From a more physical point of view, it 
can be argued that, even in a homogeneous turbulent fluid flow, particles are not necessarily 
homogeneously distributed in the whole flow domain but tend to concentrate in some flow 
regions (this is the particle preferential-concentration effect 15l-ll8|h As we are addressing 
model expressed in the framework of Reynolds-averaging approaches, the issue of whether 
this effect is specihcally well-reproduced in present formulations is not directly addressed. 
However, it is worth pointing out that particle preferential concentration effects can still 
be accounted for in the criterion (P-4). This can be achieved by considering that the same 
requirement be valid, provided that the mean value of the fluid dissipation rate is replaced 
by the local value, say (e^), ‘seen’ by particles in the specihc how regions where particle tend 
to be located. In that sense, the form of the criterion (P-4) is general enough and is relevant 
for the statistical models considered in the present analysis. 

The criterion (P-5) refers to the situation when particle inertia can be neglected (thus 
when the particle relaxation time scale becomes negligible with respect to the huid turbulence 
characteristic timescales introduced in the next section) but keeping a non-zero volumetric 
fraction Op. In that case, it is expected that particles behave as a passive (but non-vanishing) 
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scalar and the meaning of the criterion (P-5) is to reqnire that the resnlting model be 
consistent with classical scalar modeling which will be recalled in section HVBl This criterion 
complements the issne of no-spnrious drift effects, or well-mixed condition, (which will be 
shown to form one of the basis of the criteria (P-3)) and corresponds to the physical sitnation 
where particle-tracers are injected so as to create a non-homogeneous distribntion. It will be 
demonstrated in section El that this criterion is helpfnl to point to acceptable forms of the 
retnrn-to-eqnilibrinm term which typically enters the closnre of the drift vector in Eqs. (jl]). 

The criterion (P-6) is directly related to the Eqnilibrinm Eulerian Model (EEM) proposed 


for small Stokes nnmbers 
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45| (to be dehned below), which is also relevant in the present 


Lagrangian point of view. The formnlation expressed by (P-6) can be worked ont as follows. 
We consider the case of a constant particle relaxation timescale Tp. Then, starting from an 
initial condition at t = 0, the particle momentnm eqnation, Eq. fl2bjl . can be integrated to 
give 


Up(f) = Up(0)e-*/"^ + gTp (1 - / e‘'/"^U,(t') dt' 

'^p Jo 

= Vsit) + (Up(0) - U,(0)) -f gTp (1 - 


(19) 


dVs{t') 


) dt; (20) 

dt' ^ ’ 


where the second line follows from a simple integration by parts. The EEM is obtained 
by considering small Tp so that the memory of initial conditions can be neglected and 
the derivative in the integral can be approximated by its valne at time t (in other words, 
dUs{t')/dt' ~ d\Js{t)/dt since the integrand is a highly-peaked fnnction when 1). This 
yields that 


Up(t) ~U,(f) + W, 


Tr, 


dVsjt) 

dt 


( 21 ) 


where Wg = Tpg is the particle settling (or terminal) velocity. It is important to realize 
that the derivative d\Js{t)/dt is the time derivative of the flnid seen along discrete particle 
trajectories and not the flnid derivative along the flnid particle trajectory located at the same 
position as the discrete one at time t written as D\J/Dt (which is the Enlerian notation for 
the flnid particle acceleration D\J/Dt = d\J/dt + '[J ■ VU). By c omp aring with the classical 
formnlation of the EEM (see Eq. (1) in Balachandar and Eaton ISj), it is seen that this is 
eqnivalent to stating that, for small particle relaxation timescale Tp, we have 

dVs{t) DV 


dt 


Dt 


Wg-VlJ . 


( 22 ) 
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While the EEM is mostly expressed as a relation for discrete particle velocities, it is actually 
a model for the underlying velocity of the fluid seen in the limit of small Tp. At this stage, 
three remarks can be made to connect this relation to a workable criterion in our context. 
First, as we do not consider the added-mass force in the particle momentum equation, the 


parameter f3 appearing in the EEM 


18 


45| is here zero. Second, the expression of the EEM 


in Eq. is usually written in terms of non-dimensional quantities based on Kolmogorov 
scales (i.e. the Kolmogorov velocity and time scales, Up and Tp respectively), which gives 


~ - Stp 




(23) 


with = Up/ Up (the same scaling is used for and Wp""^), = f/r^, and where 

Stp = Tp/jp is the Kolmogorov-based Stokes number, which is a measure of particle inertia. 


The choice of the Kolmogorov scales is indeed relevant in DNS studies 


18 


45| . However, this 


scaling is not appropriate in the present context where we are considering stochastic models 
for high-Reynolds number turbulent flows in which a part of fluid particle acceleration is 
replaced by a white-noise term. For our purpose, it is best to introduce the timescale 
Tl oi fluid velocities (the integral timescale), which will be dehned in section IIV Al and 
to refer to the Stokes number dehned by St = Tp/T^ instead of Stp. Then, the loose 
statement of ‘small particle relaxation timescale’ can now be properly expressed as meaning 
that -C 1 or Tp -C Tl. Third, the above formulas have been obtained by considering that 
the derivative of the velocity of the huid seen is a sufficiently smooth function. This is not so 
when stochastic diffusion processes are used and when d\Js{t)/dt is white-noise. Using the 
general form of the stochastic model given in Eq. flT^ . the correct expression of the discrete 
particle velocity in the limit of small St number can be properly expressed as 

Up(f) ~U,(f) + Wp-e-*/"^ [ e*'/""D,(T) dT - [ e*'/"^B,(T) dW(T) (24) 

Jo Jo 

from which, using the same approximation, we have 

Up(f) ~ U,(f) + Wp - rpD,(f) - [ e*'/"-B,(T) dW{t') . (25) 

Jo 

If the coefficients Bg^ki of the diffusion matrix can be frozen at their value at time t, then 
the random term in this equation can be simulated (in a weak sense) as a sum of Gaus¬ 
sian random variables of zero mean and variance equal to Bl^iTp/2 [l — Yet, in 
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order to complement the criterion (P-4) which corresponds to the diffusion coefficient, we 
concentrate here on the drift vector. By comparing with Eq. ([22]), we obtain the resulting 
form of the criterion by taking the conditional average (conditioned on a given value of 
Z = (xp. Up, Us)), which gives 

(D, I Z) ~ I Z) + W, . V(U) (26) 

where we have assumed that, for high Reynolds-number turbulent flows, fluid velocity gra¬ 
dients are not strongly correlated with the (large-scale) fluid velocities so that their mean 
conditional values can be taken as being equal to the average ones. By using the Lagrangian 
formulation of DU/Dt expressed by Eqs. ([5]), we get therefore that, for St <C 1, the drift 
vector for the velocity of the fluid seen should be such that we have 

(Ds I Z) ~ (D I Z) + Wp • V(U) (27) 

where D is the drift vector used for the increments of fluid particle velocities. This is the 
basic form retained for the criterion (P-6). A weaker form can also be considered by taking 
the unconditional average, in which case we have that, for small St number, 

{D,)~(D) + Wp-V(U) . (28) 

The interest of this weaker form, which is referred to as (P-6bis) from now on, is that, for high 
Reynolds-number turbulent flows, the mean value of the drift for fluid particle velocities is 
an exact result and is equal to the fluid mean pressure-gradient (as recalled in section [TV A|) . 
Thus, the criterion (P-6bis) implies that, for small St number, we should have that 

(D,) ~ -^V(P) + Wp ■ V(U) . (29) 

Although, Eq. fl^Tj) truly embodies the criterion (P-6) in the list of requirements, both forms 
will be considered and discussed in section El 

IV. ANALYSIS OF STOCHASTIC MODELS FOR SINGLE-PHASE FLOWS 

In this section, different formulations in terms of instantaneous, fluctuating or normalized 
fluctuating velocities are considered to demonstrate the interest of the criteria selected in 
section IIII A1 
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A. Different formulations of single-phase flow stochastic models 


As Lagrangian stochastic models for fluid particle velocities attempt at reproducing some 
key (one-point) statistical properties of turbulent flows, the starting point is naturally pro¬ 
vided by the Navier-Stokes equations. Following the classical notations that use a superscript 


■*“ to indicate that we are dealing with the exact equation for a fluid particle [l| and using 
Reynolds decomposition into mean and fluctuating parts (leaving out the mean viscous term 
i>A{Ui) at high Reynolds-number numbers), the Navier-Stokes equations write 

dxt 


dt 

dU* 

dt 


= U^ 


l d{P) Y A dp 


p dxi J 


pdxi 


+ {uAui 


(30a) 

(30b) 


to model 

where the fluctuating pressure-gradient and viscous terms are to be modelec 
A reference stochastic model is the GLM (Generalised Langevin Model) 


22 | which 


represents fluid particle velocities by a stochastic diffusion process whose general form is 


dxi = Uidt (31a) 

df/j = ——- df + A df + v/ Co{e)dWi. (31b) 

P OXi 

The drift coefficient Di is usually a function of the difference between instantaneous and 
mean velocities at the particle location and is modeled as 

D. = G., (G, - X,)) = - (5 + A R(f/, - (U)) + G% (Uj - {Uj)) (32) 

where the matrix Ct- represents anisotropic effects and is subject to the condition that 
Tr{G°‘R) = 0, with Rij = (uiUj) the Reynolds-stress tensor. The drift vector can also 
be expressed so as to bring forward the timescale T^, which is a measure of the integral 
timescale of large-scale velocity fluctuations. 


A = 


Uj - m 

Tr 


+ G“(G,-([/,)) with n = 






(33) 


In these equations, the mean terms are to be understood as being the values of the cor¬ 
responding mean helds at the particle location, for instance (Ui) = {Ui){t-,x.) and, in the 
context of the present study an important point is that the GLM is formulated in terms of 
instantaneous fluid velocities. 
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So far, the GLM has been used as a reference but the present discussion is not limited 
to this model. Indeed, the GLM is one example of models written in terms of particle 
instantaneous velocities and which can be formulated as 


dxi = Ui dt 

ld{P) 


dUi = - 


p dxi 


dt + dMi{t] X, U) 


(34a) 

(34b) 


where M stands for a stochastic model expressed in terms of the state-vector Z = (x, U) 
based on instantaneous velocities and dMi its increment over a small time interval dt. The 
general model M, whose precise form is irrelevant for the present concern, is a model ac¬ 
counting for the fluctuating pressure-gradient and viscous term, as expressed in Eq. fl30|) . and 
is simply assumed to satisfy the requirements that {dMi){t, x) = 0 and {uiO dMi) / dt = —2 e. 
These two constraints translate the fact that the mean ‘force’ dMi is zero while the mean 
value of the work performed by this force is the mean dissipation sink term. In the typical 
case where dMi involves a stochastic process, the last expression has been written using 
Statonovich stochastic calculus for the sake of simplicity. In most modeling proposals, dMi 
is represented by one diffusion stochastic process which can also depend on statistics de¬ 
rived from the set of particles as well as on external helds. However, other models with 
different structures can be considered j^, 48]: for instance, dMi can represent stochastic 
models which rely on a two-level stochastic description where a random succession of elemen¬ 


tary diffusion processes are governed by a parent Poisson process 


49| . From the standard 


methodology recalled in section IIIBl it is clear that all models belonging to the class repre¬ 
sented by Eqs. (|34|l meet the criteria listed in section IIII Al and, therefore, can be regarded 
as satisfactory stochastic models for single-phase flows. 

In the passage from homogeneous turbulence to general non-homogeneous flows, other 


modeling roads have been followed. In particular, several attempts 


50 


52l | have been made at 


formulating a general model by stating that the model developed in homogeneous situations 
corresponds, in the inhomogeneous case, to the model for the fluctuating velocity, Uj, and 
that the full model for the instantaneous velocities is simply obtained by writing Ui = 
{Ui) + Ui- Other attempts |53l457l| have been made by considering that the model developed 
in homogeneous flows retains a similar form in inhomogeneous cases when it is expressed 
for the normalized fluctuating velocities, Mj/cT(j) where (T(i) is the standard deviation for 
the corresponding velocity component, and saying again that instantaneous velocities are 
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retrieved by adding the mean value at particle location (suffixes in brackets are excluded 
from the summation convention). Therefore, these formulations differ by the choice of the 
variable on which they act (instantaneous or fluctuating velocities) and by the structure of 
the model. The differences induced by these choices are now discussed. 


1. Fluctuating velocity 


To work out the relations between formulations in terms of the instantaneous velocities 
and fluctuating parts, it is convenient to start from the general form introduced in Eqs. fIMl) . 
The Lagrangian derivative is equivalent to the material derivative and for a quantity 0(t) 
sampled from a held \k(f, x) along a particle trajectory, 0(t) = T(t, x(f)), we have 


d(j){t) 

dt 


dt ^ dxk 


(35) 


Then, by writing 

Ui = Ui- (Ui) , (36) 

the evolution equation for the huctuating velocity along the same particle trajectory can be 
obtained from Eq. (I34b|l 

dui dUi d{Ui) 


dt dt dt 

dui_dU, fom ^ d{U,)\ 

dt dt \ dt ^ dxk dxk ) 

dui 

dt 


dUj ^ ld{P) ^ d{u,Uk) dm 

dt p dxi dxk dxk 


(37a) 

(37b) 

(37c) 


The hnal form is obtained by using the increments of the instantaneous huid velocities in 
Eqs. flHTl) and this shows that any stochastic model formulated in terms of U as 


dxi = Ui dt 


p dx^ 

no spurious drift 


dt + dMi{t] X, U) 


(38a) 

(38b) 


model 


is equivalent to a stochastic model in terms of u formulated as 35 


dxi = {{Ui) + Ui) dt 
d{uiUk) 


dui = 


■ dt 


dxk dx k 

(a) no spurious drift (b) production term 


Uk 1 dt +(iMj(t;x, u) 


(39a) 

(39b) 


model 
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In the formulation in Eq. fl38bp . the hrst term on the rhs (right-hand side) is the mean 
pressure-gradient and, when this mean pressure is such that the mean velocity held sat- 
ishes the divergence-free condition for incompressible hows, this ensures that an initially 
uniform (huid) particle concentration remains uniform [l, 33|. In the literature devoted to 
Lagrangian models, failure to maintain such a uniform concentration has been referred to as 
the “spurious drift ehect” or even “the well-mixed condition problem” in some works . In 
models formulated in terms of instantaneous velocities, the presence of the mean pressure- 
gradient term is evident and the issue of spurious drifts is thus trivially avoided. However, 
the situation is somewhat more involved in the formulation in terms of huctuating velocities, 
Eqs. 0391) . The hrst term on the rhs of Eq. (I39bj) . which involves the spatial derivatives of the 
Reynolds-stress tensor, must be present if spurious drifts are to be avoided. This was not the 


case in the hrst modeling attempts 
a theoretical point of view 




5l| and, although the situation was analyzed from 


33| and demonstrated from a numerical point of view [3^, this 


fact remains sometimes missed in dispersed two-phase how applications 59l-l6l|. This has 
led to a blurred vision of the rather simple issue of the spurious drift ehect and to the notion 


that hawed formulations should 
so-called drift-corrected models 


De saved by the addition of ‘corrected mean terms’, yielding 


29| . These analyses are limited to checking that the correct 


form of the Reynolds equation is satished whereas, in the context of the present study, the 
consistency of model proposals with the full structure of Reynolds-stress equations has been 
explicitly raised in the criterion (F-3). This implies that all the terms appearing on the rhs 
of Eq. (I39b[l must be addressed. 

Indeed, even if the gradient of the Reynolds stress tensor is properly introduced, the 
second term on the rhs of Eq. fl39bD is just as compulsory. This is seen by deriving the 
corresponding transport equations for the Reynolds stress components (uiUj) which reads: 


d{uiUj) , ^d{uiUj) 

~i)r~ + 

'■-V- 

{d{uiUj)) 


dxk 


S, 




dm 

dxk 


(UjUk) 


V 

correct production term 


dm 

dxk 


+ 



dm + ° 


dM,) (40) 


With the two constraints satished by dMi (see section ITV Ap . it is seen that the second term, 
labelled (b) in Eq. fl39bll . is essential to obtain the correct production term in the Reynolds- 
stress equations. Failure to account for this huctuating term means that the important 
production term in the Reynolds-stress equations is either missing or badly calculated. Such 
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a formulation would then be inconsistent with the correct form of the Rij — e equations and, 
following the set of requirements given in section IIII A1 would not be acceptable. 

Finally, once all terms are correctly handled, it is useful to compare the corresponding 
effort which is required in Eqs. fl55]) and Eqs. fl5^ . A stochastic model for U requires 3 
gradients (for V(P)) while the same stochastic model, written for u, requires 27 gradients 
(for V(U) and V{uiUj)). This represents a considerable amount of additional complexity 
and computational effort in practical calculations. For all the reasons put forward above, 
it appears therefore that formulations made in terms of the instantaneous velocity are both 
the easiest and the safest road for the construction of practical models. 


2. Normalized velocity 


Other attempts 6^ at going from the homogeneous situations to non-homogeneous ones 
were made through so-called normalized Langevin models. Recent versions of the normalised 


Langevin approach can be written as 
d 


53 


y, 


6^: 


Uj 

CT{i) 


Hi 


dt -|- 



dW, + ^ 

Tl OXk 


{UiUk} 


dt 


(41) 


Tl cr(j) V Tl oXk \ c (p 

In the following, bracketed indexes are used for the standard deviation cr(j) of the fluctuat¬ 
ing velocity component Ui to indicate that such indexes are excluded from the summation 
convention. 

Before going into the analysis of this formulation in general inhomogeneous situations, 
we hrst consider homogeneous turbulence where the model becomes 

Ui \ Ui , / 2 


d 


<T(i) 




dt -|- 



dW,. 


Tl 


(42) 


Hi 


We can derive the equations for the fluctuating velocities u* = au) - by writing 

dui = —{da(i)) + a(i)d (— 

T{i) \^ii) 


Hi 

<T(i) 


at dxk J 


(43a) 

(43b) 


Since dai^i^/dxk = 0 for homogeneous flows, this hrst proposal of the normalized Langevin 
model corresponds to the following equation for the huctuating velocities: 

1 5cr(p 1 \ r2 


dui = 


(T{i) dt 


Ui dt + au) \ — dW. 
tlJ \ Tl 


(44) 
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For stationary isotropic turbulence (where af = 2/3 fc is a constant), the normalized 
Langevin model yields that 

dui =-—dt + dWi. (45) 

'Tl y tl 

The timescale tl represents the timescale of velocity fluctuations in the stationary case and 
can be expressed as a function of k and e by introducing a constant Cq 


4 k 
3 Co (e) 


from which it results that the equation can also be written as 


(46) 


dui = - - dt + \/ Co (e) dWi 

tl 


(47) 


and since in the stationary case 2|, 1^, it is seen that the present normalized 

Langevin is identical to the SLM. Second, for the case of homogeneous isotropic but decaying 
turbulence (such as grid turbulence), the normalized Langevin model can be written 

1 daf^ 


dui = — — 


tl dt 


Ui dt “h -y/ Cq ^6) dWi 


(48) 


where the diffusion coefficient has been re-expressed with the constant Co, as just shown. 
For decaying isotropic turbulence, we have that 

1 dafi) 


2cr2) dt 


1 dk 
2 k dt 


(£) 

2k 


(49) 


and the equation for the instantaneous velocities (which, in that case, are identical to the 
fluctuating ones) becomes 


dU. = 




(50) 


By adding a possible term involving a matrix G“ as in the previous section, it is seen that, 
for the two special cases considered here of stationary and decaying isotropic turbulence, 
the normalized Langevin model retrieves the form of the GLM given in Eqs. fl3T|) - (l33|) . 

However, the situation is quite different as soon as we move out of these two simple cases. 
First of all, the formulation in Eq. fj42l) is inconsistent when it is applied to homogeneous 
anisotropic turbulence since results depend on whether the axis of the reference system 
are aligned with the principal axis of {uiUj) or not. This can be traced to the fact that 
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the model in Eq. fH2|) does not respect the criterion (F-1) (actually, the correct way to 
dehne normalized velocities should have been Ui = Uj where is the matrix 

that stands for the Reynolds stress inverse square root = 1 in matrix notation). 

Second, for general inhomogeneous flows, it appears from the form expressed in Eq. (jS]), 

lat such a formulation suffers from spurious drifts. 


, / , \ / diuiUk) , ,.11 

that {dUi) 7^ —--which shows t 


dxk 

This was recognized in earlier versions 


62| but attempts to correct the model formulation 
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57| , such as the last 


have been made mainly through the addition of ad-hoc mean terms 
term on the rhs of Eq. flTT|) . 

We can now consider the complete version given in Eq. fITT]) in non-homogeneous turbulent 
flows. By following the same derivations as in Eqs. fH3|l . this normalized Langevin model is 
equivalent to the following equation for the fluctuating velocities 


diUiUk) 1, j 1 1 ^cr(i) , jj 1 , 

dUi = —;;- dt —{UiUk} -— dt + UiUk -— dt ■ 


dxk 


cr(i) dXk 


■(Tp) dxk 


.J_?^u,dt-^dt + a,.f^dW 

m tl tl 


^‘‘production term” 


Langevin term 


(61) 

diu'U ^ 

Since the constraint that (dui) = ——- dt is now enforced, there is no spurious drifts 

OX]^ 

and the model is consistent with the Reynolds equation, indicating that the hrst criterion 
put forward in section UlI Al is met. However, as seen by comparing Eq. (I5T|) with the exact 


equation for fluctuating velocities in Eq. fl39bp . it is evident that the exact term —Uk 


dm 

dxi 


IS 


not retrieved. This means that the essential production terms in the corresponding Reynolds- 
stress equations are mishandled and, therefore, that this model is inconsistent with the 
correct structure of — e equations, as presented in Eqs. flTTl) . In other words, present 
normalized Langevin models do not satisfy the criteria (F-3) and the failure to reproduce 
the correct structure of Reynolds-stress equations is a severe limitation. 

It is also worth noting that, if the previous equations have been developed using the 
simplest form of the Langevin model, these results does not depend on the special form of 
the normalized model. Indeed, if we consider a similar model written as 

d ({uiUk)' 


a(i)J dxk 


(X(i) 


dt + dMi{t, X, u) 


(52) 


where dMi is expressed as a function of the normalized fluctuating velocities u which, in 
present formulations are (wrongly) dehned as Ui = Ui/a{i), and is such that {dMi) = 0 
and {ui dMi) = 0 to respect stationary isotropic turbulence conditions, we have for the 
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fluctuating velocity an equation similar to Eq. fl5T|) . namely 

, dtuiUh) , , ,1 da(i\ , 1 da(i\ , 1 da(i\ , , 

dui = — dt — (uiUk) - —- dt + UiUk - —- dt H- Tr^Ui dt + a^t) dMi. (53) 

^ il) 


dxh 


o-(i) dxk 


'(T(j) dxk 


Thus, whatever the chosen form of the stochastic model dM^ and regardless of variants in the 


expression of the diffusion matrix 
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3, 


the same conclusions about the lack of consistency 


and the failure to respect the criteria (F-3) still hold, making these models unacceptable 
descriptions of turbulent flows. 

B. Consistency with scalar modeling 


Similar issues exist for the formulation of scalar models and, furthermore, consistency 
with classical scalar modeling is explicitly used in the criteria (P-5) retained in the list of 
requirements set forth in section UlIBI for two-phase flow models. It is thus worth discussing 
the relevant points of classical scalar modeling. 

For this purpose, we consider a passive scalar 0(f, x) transported by the flow and which 
is the solution of the exact advection-diffusion equation 


d(j) d(j) _ 
dt ^ dxk dxl 


(54) 


where P is the scalar diffusivity. Assuming, still for the sake of simplicity, that we are 
dealing with high Peclet-number flows, the mean diffusion term can be neglected and the 
mean scalar equation (using 0 = (0) -|- 0') has the form 


= 0 


( 66 ) 


dt dxk dxk 

where {uk 00 the scalar flux which must be modeled. In second-order turbulence modeling, 
non-local closures for the scalar flux are obtained by considering the transport equation for 
the scalar fluxes which write 

^M') , , d{UiUk(l)') ! ^0(0) / , ,T, 

~dr + dx, “ 

where Tj represents the correlation (0' + A0') which is usually modeled 

as a function of the scalar fluxes, for example T* = —Bikiukcf)')■ 

A classical approach is to consider the algebraic relations which result from the second- 
order transport equations when convective terms are neglected 6^. This approach was 
























developed to show that, in that case, one should obtain expressions consistent with the 
general diffusivity concept. Indeed, neglecting transport terms in the above equation for the 
scalar fluxes gives 



( 57 ) 


which yields that 




djcj)) 

dxj 


(58) 


d{Ui) 


with Rij = (uiUj) and Oij = Bij + . 

This behavior has been found to be realistic both in numerical and experimental 
investigations 64l. l65| . It is thus important to check that various formulations can yield 


similar expressions when the same hypothesis are made. For example, starting with the 
formulation in terms of instantaneous variables that was already used at the beginning of 
section IIV Al this consists in adding the instantaneous scalar value attached to a particle 
trajectory in the state-vector Z = (x, U, 0) and in replacing the exact equations by modeled 
ones. Using for instance the GLM already introduced in section IIV Al for the dynamical 
variables and adding a model for 0, we obtain a typical model structure as 


dxi = Uj dt 


dUi = dt + Didt+ JCo{e)dWi 

P OXi 

d(j) = Afj, dt. 


( 59 ) 

( 60 ) 
( 61 ) 


In the last equation, stands for a model for F A0' and is subject to 
{A^ I U) = 0. This is referred to as the micro-mixing modeling issue j^, 3 


re constraint that 
5|] which has been 


the subject of ongoing research efforts. For the sake of our present discussion, it is sufficien 


to use the simple lECM (Interaction by Exchange with the Conditional Mean) model 
which reads 

0-(0|U) 


^1 


A^ — 




(62) 


where is the scalar mixing time scale. Then, using either a proper PDF derivation pij or 
a short-cut method as in section IIV Al it is straightforward to show that the scalar fluxes 
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are the solutions of the transport equations 


d{ui 0 ) ' d{ui 0 ) d{ui Uk 0 ) 

dt ^ dxk dxk 


(UiUk) 


dxk 


{Uk 0 ) 


dm 

dxk 





{Uk4>) (63) 


In this equation, the last term on the rhs represents a model for dtj resulting from the 
specihc choice of the GLM for particle velocities and the lECM model for 0. Other closures 
would result in different expressions but with a similar form. The important point in the 
present context is that the formulation in terms of instantaneous variables gives the correct 
structure and is, thus, quite consistent with the asymptotic analysis developed above to 
derive algebric relations when transport terms are neglected. 

Should alternative formulations be considered, for example in terms of fluctuating scalar, 
the developments presented in section IIV Al indicate that the evolution equation for 0' is 


d0' 


OXk 


(9(0) , . , 

-Uk^-^dt +A^>dt 

0Xh 


(a) no spurious fluxes (6) production term 


(64) 


where corresponds to the micro-mixing model retained and expressed now in terms of 

the fluctuating components. For instance, for the lECM this fluctuating term is given by 

A^i = A ((0 I U) — (0))/r(^. The hrst term, labeled (a) in Eq. (IMD . means that 

the correct form of the mean scalar equation, Eq. flS^ . is retrieved and is necessary to 

avoid spurious scalar fluxes. The second term, labeled (6) in Eq. fIMl) leads to the correct 

production term, — (ujUfc)——, in the transport equation for the scalar fluxes in Eq. (I6^ 

OX]^ 

and, without this term, the asymptotic analysis that retrieves the scalar turbulent diffusivity, 
as in Eq. fl58|) . breaks down. 


V. ANALYSIS OF STOCHASTIC MODELS FOR TWO-PHASE FLOWS 

In this section, we turn to stochastic models used for the simulation of fluid velocities 
seen by discrete particles. The discussion is developed by considering some specihc models 
which have been proposed in the literature. However, it is worth repeating that the aim is 
to help clarifying how formulations can be assessed, what issues are misleadingly taken as 
important and which ones should be carefully addressed when building new ideas. 
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A. Different formulations of two-phase flow stochastic models 


1. First Langevin proposals 


One of the first proposals to introduce a Langevin model for the velocity of the fluid seen 
relied on a formulation based on the instantaneous velocity along particle trajectory 19 |. 


We discuss this first proposition since it is interesting to illustrate the stochastic modeling 
procedure, in spite of some limitations outlined below. The equation was written in a discrete 
form but, using present notations for the sake of consistency with the other formulations, 
this model for U^, can easily be expressed as 

dUs,i = dt + iyA{Ui)dt+ (f/p,fc - Us,k) dt + G*, {Us,k - {Uk)) dt + B dW, (65) 


p dxi 


dxi 


In this equation, the matrix is defined by the two timescales T£|| and T£_|_, which 
correspond to Csanady’s expressions (though the precise form of these Csanady’s formulas 
were not given in the original work) for the timescales of the fluid seen in the directions 
parallel and perpendicular to the direction of the mean slip velocity (Ur)m = (Up)m — (U) 
respectively, and is given by [19 1 


G * _ 

ik — 


n~'* 

-^L,_L 


dik 


1 


T* 


T* 


Ti Tk 


( 66 ) 


■L,\\ ^L,±_ 

with Ti = Ur^i/\ {Ur)m \ the normalized vector aligned with the mean slip velocity. The 
Csanady’s formulas for the timescales || and T£ are jll 

Tl Tk 


Tlw = 


'l + /d 


2l(U.), 


TU = 


2k/?, 


'1 + 4/3 


2l(U.). 


(67) 


2k/? 


where the expression of is given in Eq. 


and where j3 = T^/Te is the ratio of the 


Lagrangian integral time scale to the Eulerian one Te- 

For two-phase flow modeling, this article was interesting in that it was one of the 
first attempts to build bridges between the single- and two-phase flow situations. Another 
interest is that it introduced the notion of the consistency between the resulting model for the 
drift velocity and classical scalar modeling as detailed in section HVBl This consistency 
issue has been retained in the present list as criteria (P-5). It is seen that in Eq. fl65l) the 
second term on the rhs of the equation stands for the mean viscous term. In the limit of 
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vanishing particle inertia, this model gives an expression similar to the SLM to which the 
mean viscous term has been added. Such a formulation is indeed consistent with the low 
Reynolds-number form of the Reynolds equation but it must be noted that it would not give 


the correct low Reynolds-number form of the second-order equations 


42 


43| . However, 


m 


the present context, we have limited ourselves to high Reynolds-number flows and, therefore, 
the mean viscous term can be neglected. 


This hrst proposal has often been cited as a ‘stochastic model’ 
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68| to be com¬ 


pared on an equal footing with the more complete formulations that will be discussed below. 
However, for all the interests of the form expressed by Eq. ([H5D, this is a misleading pre¬ 
sentation and an incorrect statement with regard to its actual status. Indeed, the diffusion 
coefficient, written as B in Eq. fl65|) . is never specihed 19|]! In the context of the orig inal 
study, this is not surprising since, as it transpires from the very title of the article [l9|, its 
purpose was to derive a set of mean-held equations from this stochastic description. As the 
set of mean-held equations was chosen to be limited to the particle mean velocity and the 
huid-particle correlations, the expression of the dihusion coefficient is not required. If such 
a proposal leads to as a realizable Eulerian model, it must be stressed that, in no way, can 
it be called a Lagrangian stochastic model. From the list of criteria in in section IHI Bl it 
is indeed obvious that the criterion (P-2) is not met. In the particle-tracer limit, the form 
given in Eq. fl6^ reverts to the form of the SLM and, in that sense, all the terms appearing 
on the Ihs of Eqs. mm are obtained. However, since the dihusion coefficient is unknown, 
it cannot be assessed whether the corresponding terms Tjj on the rhs of Eq. (II Teh is such 
that l/2Tjj = —e with the constraint that e > 0. Therefore, the criterion (P-3) cannot be 
checked. This is also true for the criterion (P-4) which is directly related to the dihusion 
coefficient and is therefore not met, whereas the criterion (P-5) is satished. In order to 
assess criterion (P-6), we consider the limit of non-vanishing but small Stokes numbers. In 
that case, we can assume that the mean relative velocity {Ur)m is equal to the particle 
settling velocity introduced in section IIIIBl that is (Ur)m — W^. Then, given the 
expression of the Csanady’s timescales in Eqs. (1671) . which can be re-written under the form 
T£ II = Ti/\/l + with X = St X {g T^)/y/2/3k, we have to hrst order in St that || — 
as well as T£ ~ T^. This means that the mean conditional average of the drift vector 
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(Ds I Z) which appears in the criterion (P-6) is such that 


{Ds,^ I Z) 


1 d{p) Us,i - m d{u,) 

p dxi Tl dxk 


(68) 


The hrst two terms on the rhs correspond to the mean conditional increments for fluid 
particle velocities modeled with the SLM. We can thus re-express this equation as 


I Z) 


{Di I Z) -|- Ur,k 


dm 

dxk 


(69) 


where D is an acceptable model for the drift vector for fluid velocities, as demonstrated 
by the analysis carried out in section IIV A[ By comparing with the formulation of the 
criterion (P-6) in Eq. fl27D . we can see that the second term is not exactly retrieved since 
the instantaneous relative velocity appears instead of the mean one. Thus, strictly speaking, 
the criterion (P-6) is not satished, though it is also readily seen that the weak form (P-6bis) 
is respected. 

In summary, the proposal in Eq. (I65|) must be regarded as an incomplete Langevin model 
and cannot be retained as a proper PDF description. At this stage, it is an often-seen temp¬ 
tation to close the diffusion coefficient by retaining the closure used in the GLM approach for 
single-phase flows, that is B = Cq (e). Yet, it can be seen from Eq. (|65l) that the return- 
to-equilibrium term is not based on a scalar timescale (as in the SLM with T^) but has a 
non-isotropic form. In other words, using the fluid-limit value of the diffusion coefficient 
would clearly lead to a violation of the criterion (P-4). This central point in the formulation 
of stochastic models for the two-phase flow situation was only recognized afterward j4|. 


2. Complete Langevin models 


The velocity along particle trajectory was later called the velocity of the fluid seen and 
models were written as proper stochastic differential equations. The general PDF framework 
was set up jd] and a new Langevin model was developed as one example of the general 
methodology. For the sake of a simpler presentation, we assume here that the mean slip 
velocity (Ur)m is aligned with the hrst coordinate axis and, based on this assumption, the 


33 











first complete Langevin model is given by 


dt + {{Up,u)m - {Us,k)m) ^ dt 


P dxi 


U s^i 

rji^ 


dxk 


dt+ \l (e) ( Cobi^i)kjk + -(6p)fc/fc - 1) ) dWi (70) 


In this equation, A; is a new kinetic energy which stands for the fluid normal kinetic energies 
weighted by the Csanady’s factors and is expressed by: 


k = 


SEliMfil 

2 


h = 


Tr 


T, 


(71) 


L,i 


With the assumption on the coordinate system, the time scales T£ are equal to 

^L,l - ^L,\\ ; 


rynH' 

-‘■L,2 = -‘■1,3 = -‘l,± 


(72) 


where T£ || and T£ are given in Eqs. (ETj). The form given in Eq. (1701) in the reference system 
aligned with the mean relative velocity is best to bring out the characteristic features of the 
model and the forms of the drift and diffusion coefficients. Yet, in a general coordinate 
system, this form of the complete Langevin model is expressed by 

dt + {{Up,k)m - {Us,k)m) ^ dt + G*, {Us,k - {Us,k)J dt + B,, dW, (73) 

and where the explicit expression of the diffusion matrix 


p dx 

where is the same as in Eq. 


Bij is detailed elsewhere 




It is seen that the closure of the drift vector in Eq. (1701) differs slightly from the one chosen 
in Eq. (|65ll . However, the key point is that the diffusion coefficient (actually, the diffusion 
matrix) is explicitly formulated. Therefore, the hrst part of the criterion (P-2) is now met, 
which explains why this model has been referred to as a ‘complete Langevin model’. It is 
also clear that such a model allows practical numerical predictions to be carried out ^ and, 
consequently, the criterion (P-2) is fully satished. In that sense, this Langevin model can 
be truly referred to as a PDF description for two-phase flows. It is also straightforward to 
show that, when particle inertia becomes negligible, the model given in Eq. flTOl) reduces 
exactly to the SLM for fluid particles which is one of the models that meet the lists of 
requirements set for single-phase flows. Thus, the criterion (P-3) is satished. The specihc 
closure of the diffusion matrix is such that, for a homogeneous turbulent huid how laden 
with discrete particles, the kinetic energy of the huid seen follows the same statistical law as 
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for the fluid kinetic energy, d{l/2 (u^)m)/dt = —(e). As discussed in section ITlIBl when the 
criterion (P-4) was introduced, this is indeed what can be expected from a physical point of 
view. This apparently-trivial constraint is actually an important issue in the construction 
of complete stochastic models for two-phase flows and was stressed accordingly For our 
present concern, this means that the criterion (P-4) is also met. However, in the form chosen 
for the drift vector in Eq. d70|) . it is seen that the return-to-equilibrium term was written as 
a return of the instantaneous fluid velocity seen to its mean value at the same location. This 
results in a zero-contribution term in the corresponding mean equation for the drift velocity 
and, therefore, the criterion (P-5) is not respected by the model given in Eq. flTO]) 691. 

This limitation was later overcome by expressing the return-to-equilibrium term as a 
return of the instantaneous value of the fluid velocity seen to the local mean fluid velocity. 
This led to the hnal formulation of this complete Langevin model 23| which, using the same 


simplifying assumption as in Eq. d70|) that (U^) is aligned with the hrst coordinate axis, is 
expressed by 

dUs, = dt + {{Up,k)m - {Uk)) ^ dt 


p dxi 




m 


-L,{i) 


dxk 


dt+J{e){ C^\i^k/k + -{\i)klk - 1) dWi . (74) 


whereas the complete formula in a general coordinate system is 

dU. k = dt + {{Up,k)m - (f/fc)) dt + G*k {Us,k - (Uk)) dt + B,, dW, (75) 


P dxi 


with the explicit formulation of the diffusion matrix detailed elsewhere 231. 

Like its preceding version, this formulation fulhlls the criteria (P-1) to (P-4). However, 
with the new form for the return-to-equilibrium term in the drift vector, the corresponding 
equation for the drift velocity is now 

1 


“T \^p,fc/r 


d \ d 

[af{UiUk)] — y— (”^6) 


ttp dxk 


r\, I X^U.K/rn rN 

at oxk afOXk 

where and ap are the fluid and particle volumetric fractions (with aj + ap = 1). In 
Eq. H76]l . the fluctuations of Ug^i are dehned as uG = Ug^i — {Ui) (Note that a specihc 
notation is used for the ‘fluctuation’ of the velocity of the fluid seen while usual notations 
are kept for particle velocities and for the fluid velocities within the fluid phase, see the 
discussion on the dehnition of fluctuating velocities in section IV A 3^ . It is important to 
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stress that, in the scalar limit case (when the particle relaxation time Tp is negligible with 

respect to the fluid turbulence time scale Tl but with a non-vanishing volumetric fraction), 

Vd^i reduces to the ordinary turbulent correlation between the particle concentration and the 

fluid velocity, which is the case discussed in section IIVBI By making the same assumptions 

as in section IIVBI this equation yields that 

_{ I d Id 1 

= {G (77) 

L ^ f ^ ) 

which is in line with classical scalar modeling as recalled in section IIVBI Indeed, for the 
case of constant correlations, we have that 


f 1 1 ^ 

V,, = -Kknp,i)m + -{UkUl) ^ (78) 

CXp Oif j OXi 

which is similar to the result given in Eq. flSSp . showing that the criterion (P-5) is fulfilled. 
For the analysis of the criterion (P-6), we can draw on the developments obtained in sec¬ 
tion IV A II in the limit of small Stokes numbers. Thus, to first order in St, we have that 
T£|| ~ Tl, T£ _|_ ~ T^, as well as ~ 1. Using the general expression of the drift vector 
in the complete model in Eq. fl75|) . we obtain now that the mean conditional average of the 
drift vector (D^ | Z) is 


,,, iw u.,-m , , am 

\^s,i I --^ mr,k)7 


(79) 


O m ' X'-'T,K m r\ 

P OXi Tl dXk 

Using again the drift vector of the SLM for fluid particle velocities and the fact that, for 
St 1, (Ur)m — Wg, this is equivalent to 

dm 


{Ds ,7 I Z) = (A I z) + Wp,k 


dxu 


(80) 


This is indeed the form stated for this criterion, as expressed in Eq. ([27]), showing that (P-6) 
is satisfied, from which (P-6bis) is obviously also met. 

Consequently, the complete list of criteria of section IIIIBI is satisfied by the Langevin 
model given in Eq. flTlD or in Eq. fl75]) and, as such, this model can be assessed as being a 
satisfactory description of two-phase flows. 


3. Fluctuating or instantaneous fluid velocity seen 

So far, the models considered in sections IV A II and IV A 21 for two-phase flows have been 
presented in terms of the instantaneous fluid velocity seen. Formulations in terms of the 
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fluctuating component of this fluid velocity seen are, of course, possible but reveal themselves 
to be not only cumbersome but also trickier than in the single-phase case. One reason is 
rooted in the fact that ‘fluctuations’ must be carefully dehned since averages should be hrst 
properly dehned as meaningful quantities over the respective huid and particle phases. For 
example, two ‘huctuations’ can be dehned from the instantaneous value of the huid velocity 
seen: as in section IV A 21 one can introduce huctuations as the diherence with the local mean 
huid velocity (properly dehned as an average over the huid phase only), uO = Ug^i — {Ui), 
whose mean value (over the particle phase) is not zero but the drift velocity or the 
‘real’ huctuation, Ug^i = Ug^i — {Ug^i)m, whose mean value is indeed zero (obviously, these 
two huctuations are related through Ug^t = ^ — Vd^i. Note that {u'g^Up^k)m =J^ajUp,k)m in 

Eqs. fl75D - fl75D ). These distinctions were made clear in Eulerian descriptions jl9l. l44| and, in 
the complete PDF approach to two-phase hows, this means that the complete theoretical 


framework 


m, 


2l| outlined in section [TTl must be carefully followed. 


To illustrate this point, if we consider a stochastic formulation for the instantaneous huid 
velocity seen which has the following form: 

ld{P) 


dUg i 


p dxi 


dt + Aidt -\- G*j^ (Ug^k ~ {Uk)) dt Bik dWk 


(81) 


then, the corresponding equation for the huctuating component u'g ^ as dehned above is 


du' i = 


1 d[af{uiUk)] 


dt — [Up^k ~ {Uk)] ^ dt + Aidt + G*j^u'g kdt + BikdWk (82) 


«/ OXk OXk 

Handling directly this form of a model raises questions as to the exact formulation of the 
hrst term on the rhs of Eq. fl5^ . Indeed, an alternative formulation can be expressed as 


duU = dt - [Up^k - {Uk)] ^ dt + A, dt + dt + B^kdWk . (83) 

This dual form of the equation for the huctuating part of the velocity of the huid seen is not 
related to a particular choice of the drift and dihusion terms. For example, in recent arti¬ 


cles 


£ 3 , 


68j, it is said that for the value of the term Ai in the drift vector which corresponds 


to Eq. fl65l) . the equation for this huctuating component is 

whereas the corresponding hrst form given above would be 


(84) 


(85) 
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The difference between Eq. fl82D and Eq. fl83D (or between Eq. flHTl) and Eq. fl85|) i is related 
to whether two-way coupling effects are accounted or not. The form in Eq. flH^ (and in 
Eq. 085 p i is obtained when the mean fluid velocity is given by 

ld{P) 1 d[af{uiUk)] 


am ^ ^^^^am 


dt 


dxi 


p dxi 




dxi 


( 86 ) 


which contains the correct expression of the Reynolds stress for the fluid momentum equation 
in the two-phase flow situation { 4 , 21|. On the other hand, the form in Eq. (and in 
Eq. flHTp i is obtained by considering that the equation for the mean fluid velocity follows 

ld{P) d[{uiUk)] 


am ^ ^^^^am 


(87) 


dt dxk p dxi dxk 

Eq. 0571) is, of course, the Reynolds equation for the fluid phase treated as an incompressible 
single-phase turbulent flow (in the high Reynolds-number limit). On the other hand, in 
Eq. 086p . it is seen that a two-way coupling effect is present since the reduced volumetric 
fraction occupied by the fluid is accounted for through af{t, x). Note, however, that Eq. 086P 
represents only a partial account of two-way coupling in the sense that only the volume effect 
is included. In order to properly account for two-way coupling, momentum exchange terms 
should be added for the stochastic model for the velocity of the fluid seen and into the 
Reynolds equation for the fluid phase. To the authors’ knowledge, this has been prop erly 
proposed only for the complete Langevin model (see detailed presentations jd, 6, 
considering that the term 


m by 


A. 


P^f,i 


Pp ( ^P,i 


afp 


Tr, 


( 88 ) 


which represents the effect of the particle phase on the fluid velocity seen is added to the 
rhs of Eq. (17^ . In that case, volumetric and momentum two-way coupling are taken into 
account and the corresponding equation for the fluctuating part of the velocity of the fluid 
seen is 

dm 


duu=-mmdt 

’ af dxk 


[Up,k - m] 


dxi 


■ dt ylj dt f. dt Pj^kdWk 


Oip Pp 

afp 


U.A 


Up,i 


Tp 


ik^s,k 


Up A 


Tp 


dt . (89) 


Note that, in the usual case of a polydisperse two-phase flow, the last two terms on the rhs 
of Eq. (1591) cannot be expressed only in terms of the mean and fluctuating fluid and particle 
velocities (which is a further argument to suggest that formulations in terms of fluctuations 
are more cumbersome than the ones in terms of instantaneous values). 
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Though the purpose of our analysis is not directly about two-way coupling effects, the 
above discussion is interesting to reveal that, for a given stochastic model, the form of the 
corresponding equation for the fluctuating velocity components is not immediately deter¬ 
mined. Indeed, if we consider for example the complete Langevin model, its formulation 
depends on whether one wants: (a) to disregard all two-way coupling effects (both volumet¬ 
ric and momentum), which leads to Eq. (18^ : (b) to account for the volumetric effect for 
the fluid phase, which leads to Eq. fl82|) : (c) to include two-way coupling effects, in which 
case the model for the instantaneous fluid velocity seen, Eq. fl7iD -l-Eq. (188|) . corresponds to 
Eq. flS^ for the fluctuating component of the velocity of the fluid seen. 

Going back to the difference between Eq. (|8^ and Eq. (|8^ (or between Eq. (|Mll and 
Eq. fl85|) ). this means that an additional term, equal to {uiUk)l/cxfVap and proportional to 
the gradients of the volumetric fractions Vop, is introduced in Eq. (j83]). Depending on the 
choice made for two-way coupling effects (see (a) or (b) above), this (potentially) spurious 
term may not be negligible as particle distribution is rarely homogeneous. Consequently, 
numerical simulations which are performed from the stochastic differential equations for the 
fluctuating parts of the velocity of the fluid seen can be potentially flawed. 

Clearly, formulations in terms of so-called ‘fluctuations’ can easily become intricate with 


many terms, making their manipulation slippery w 
physical picture. Thus, contrary to some statements 


rile they do not necessarily clarify the 




68 [ |. a more practical way to account 


for the velocity of the fluid seen is to retain formulations in terms of the instantaneous 
velocity. 

4- Normalized Langevin models 


Another category of models consists in proposing Langevin-type of models developed in 
terms of the normalized fluctuating fluid velocity seen. To give one example, this normalized 


fluctua 


as [53|, 


in g ve locity (using the notations introduced in section IV A 31) is sometimes modeled 
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d( — 1 = 

V(i). 





■ dt + i — dWi + 

Tl(^({) V Tl OXk 


d fiKK) 


dt 


(90) 


CTp) / \1 + St^ 

Such formulations have been mostly used for boundary-layer simulations often in relation 


with particle deposition issues 


54l |. It is seen that these propositions include terms 
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written as functions of the particle Stokes number, which with present notation is St = 

Such models are clearly able to predict particle statistics and, in that sense, they satisfy the 
criterion (P-2). However, when the Stokes number goes to zero (that is when particle inertia 
becomes negligible), these models revert to the normalized models which have been discussed 
in section IIV A 11 As analyzed there, these normalized Langevin models do not have the 
correct transformation and invariance properties and are inconsistent with the structure of 
Reynolds-stress equations: they are, therefore, unsatisfactory descriptions for single-phase 
turbulent flows. This means that the criteria (P-3) and (P-1), as well as the two forms of 
the criterion (P-6), are not respected. Consequently, these models cannot be regarded as 
acceptable descriptions for two-phase flows. 


5. Hybrid DNS-stochastic approach 


As already mentioned above, a recent proposal introduced a new formulation for the ve¬ 
locity of the fluid seen 68|. Using the present notations, this proposal consists in simulating 
as the solution of the following stochastic differential equation 

dUg^i =-^^ dt -\- V /S.{Ui) dt {Up^k — Us,k) ~ dt A[ dt 


p dxi 


dxk 


+ G:,{Us^k-{Uk))dt + BikdWk (91) 


where a new term A'- is added to the drift vector in Eq. 0651) and is given by 


A' = 


9{'^s,i'^P,k)m diUiUk) 


(92) 


dxk dxk 

This model is expressed in terms of the fluctuating components of the velocity of the fluid 
seen and the discussions in section IV A 31 are thus relevant. In particular, the formulation 
put forward [68| is along the one expressed in Eq. (18T|) to which the extra term A' is added, 
which results in 


dUc 1 Q dt I 0 j ^s,k dt ^ikdWf^ . 




(93) 


In this formulation, the matrix G*j^ is expressed as the inverse of the matrix Tik which is the 
matrix of the decorrelation time scales of the fluid seen and given by 

dm 


^ik — m + 


ik 


dxi 


m*k = - ir)ii ,%k= {umx{<imum))rndt (m) 
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where the values of the matrix Tki are obtained as statistics extracted from a DNS of the 
same particle-laden flow that is considered. Then, the diffusion matrix is determined by the 
following equality [68|: 


- BuBki — 


(95) 


Note that the notation which was used 


68 | should in fact refer to This proposal 


introduces new elements and, since DNS results are built-in, it is referred to as an ‘hybrid 
DNS-stochastic’ approach. Some more comments are in order about the role of DNS in this 
hybrid approach. In this formulation, the matrices Tik and G*^. are not obtained hrst on 
some simple test cases (thus from some sample DNS) and then used as a predictive model 
in any geometry. Actually, for each geometry and for each flow considered, a DNS must be 
run beforehand in order to obtain the correlation matrix which is used to extract the matrix 
G*i^ fed into the governing equation for the velocity of the fluid seen, Eq. fl9T]) . This means 
that a concurrent run of a DNS and of the PDF model must be carried out. 

Compared to the propositions addressed in the previous subsections, the additional term 
A' in the drift vector is also a new element. This term is motivated by an analysis of 
the form of the exact equation satished by the drift velocity in the limit of high-inertia 
particles {St S> 1 or Tp —)■ oo), while the two expressions that yield the matrix as well 
as the diffusion matrix Bik in Eqs. flOTl) and fl9^ respectively are based on a so-called ‘local 
homogeneity assumption’. In Eq. flOTH . it is seen that the mean viscous term is present as in 
the proposal in Eq. fIBS]) . However, as already indicated in section IV A 11 this is inconsistent 
with the proper low Reynolds-number form of the Reynolds-stress equations and should be 
avoided, though it seems to be systematically presented [291]. Since we are only considering 
high Reynolds-number flows, this term is disregarded in the present analysis. 

From a physical point of view, some of these assumptions raise questions. The local 
homogeneity assumption means that the underlying fluid turbulence is assumed to remain 
homogeneous and stationary so that the turbulent characteristics seen by particles (which 
are functions of the elapsed time) can be taken as stationary processes during a ‘long-enough’ 
period (at least of the order of the particle relaxation timescale Tp) for these statistics to 
reach their local ‘equilibrium values’. For high-inertia particles which can typically cover 
large distances during a time of the order of their relaxation timescale, this amounts to 
assuming that the fluid turbulence remains homogeneous and stationary over considerable 
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length and time scales, for the hypothesis to apply. This 


similarities with one development proposed by Pope 


lybrid-DNS approach bears some 


71| . However, that study [7l| was 


for homogeneous shear fluid flows where, once correctly rescaled, fluid particle velocities 
constitute truly statistically-stationary processes. Furthermore, its purpose was mainly to 
demonstrate the potential of properly-dehned Langevin models, suggesting that there is still 
room for considerable improvement in the development of particle stochastic models j^ . 
rather than as a general methodology to apply in any flows. On the other hand, the issue 
of whether consistency results are to be expected in the high-inertia limit is worth noting 
and will be taken up in the next subsection. In the context of the present analysis, we 
now concentrate on how this new proposal stands with respect to the criteria retained in 
section IIIIBl 


Since the drift vector and diffusion matrix are explicit in Eqs. finT]) - fl^^ . the stochastic 
model is complete. From the choice of the drift term, it is also seen that this proposal 
is such that the criterion (P-5) is satished and that the fluid limit yields a model whose 
form is acceptable. In the limit of small St, the additional term in Eq. fl9^ does not 
contribute to the hrst-order development in St and, in that limit, the form of the drift vector 
is identical to the one given in the incomplete Langevin proposal in Eq. Consequently, 
the analysis performed in section IV A II applies here, showing that (P-6) is not satihed 
but that (P-6bis) is. However, the main characteristic of this proposal is that the matrix 
entering the return-to-equilibrium term in the drift vector is provided by the solution of a 
DNS for the same flow which is modeled by the proposal in Eq fl?I|) . As such, it is clear 
that the criterion (P-2) is violated and that this proposal cannot be referred to as a ‘model’. 
Indeed, as indicated above, the formulation requires that a DNS solution be performed 
beforehand for each flow to feed the drift term of a model which is supposed to provide 
predictions on this very flow. If such formulations can indeed be considered as making 
interesting consistency checks for specihc closure propositions, it cannot be accepted as a 
proper stochastic model for two-phase flow simulations. Furthermore, it appears that the 
present formulation suffers from some inconsistencies that limit its applicability. Indeed, 
if we consider one of the simplest situations, namely of homogeneous isotropic decaying 
turbulence, the present hybrid DNS-stochastic approach would consist in running hrst a 
DNS to obtain the timescales necessary to dehne the matrix in Eq. fl9T)) . from which the 
diffusion matrix is derived through Eq. fl95|) . The latter equations is the expression of the 
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classical fluctuation-dissipation theorem for stationary processes [30, 1 m 1, which means that 
we have d{ul)m/dt = 0 instead of the natural decaying law for the turbulent kinetic energy. 
Thus, the criterion (P-4) is also violated. It is also evident that the proper fluid limit is 
not retrieved as the decay law for single-phase turbulent homogeneous flows is not correctly 
retrieved, showing that (P-3) is also violated. For an approach based on DNS results, this 
is a severe shortcoming. Furthermore, it can be noted that the procedure of coup ling DNS 


and Langevin models raises consistency questions, especially in the fluid limit 


46|. 


B. Discussion on additional criteria for two-phase flow modeling 


In the list of criteria set forth in section IIIIB1 it is seen that the fluid, or particle- 
tracer {St 1), limit has been particularly emphasized (in the criterion (P-3)). It can 
be wondered whether similar consistency limits are to be selected in the other limit, when 
particle inertia is very high {St 3> 1). However, three remarks can be made. First, it must 
be remembered that the present form of Langevin models (i.e. modeling the velocity of the 
fluid seen by a diffusion stochastic process) has less justihcations than in the fluid case. For 
high-inertia particles, the well-known frozen-turbulence hypothesis can be applied and it is 
easy to show that the increments of the velocity of the fluid seen should then be governed by 
spatial correlations. In that case, ((dUs)^) scales as ((e) |Ur| over a time interval At, 

provided that the relative mean distance covered by particles Ar = |Ur|At remains in the 
inertial range length Ar L. It remains to be seen whether such limits are encountered 
as they imply the existence of a regime where |Ur| is high enough so that Lagrangian fluid 
correlations can be neglected while we still have that Ar = |Ur|At -C L. If present, such 
limits are indeed not reproduced by present Langevin formulations as discussed at length in 
the construction of the complete Langevin model j4|. This suggests that, in the high-inertia 
limit, the issue would be to devise altogether new stochastic model formulations rather than 
imposing requirements on the present ones. On the other hand, when dealing with high- 
inertia particles, a second remark is that the detailed form of a Langevin-type of model 
is less an issue if macroscopic particle statistical properties, such as diffusion coefficients, 
kinetic energies, etc., are sought. For these properties, the important element is mainly to 
retrieve the correct limit of the integral time scale of the velocity of the fluid seen in this 
limit. This is indeed what expressions such as Csanady’s formulas in Eqs. fl67)l are doing 
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without having to change the form of the Langevin equation itself. Finally, with respect to 
the particle relaxation timescale Tp and since S'f S> 1 means that T£ Tp, the fluid velocity 
seen becomes a fast-variable and can be safely taken as acting as a white-noise on discrete 
particle trajectories. In that case, the form of the diffusion matrix in Langevin formulations 
is more important than details of the drift vector (apart from the return-to-equilibrium term 
and the integral time scale). Among other asymptotic cases, a comprehensive discussion of 
this high-inertia limit was proposed jb] along with the corresponding consequences on the 
development of suitable numerical schemes 6|. 

Thus, the development of new stochastic models that would capture both Lagrangian 
statistics (time-spectrum) as well as Eulerian ones (space-spectrum) is still an open issue. 
Yet, this question must be weighted against what would be gained from such developments. 


C. Open issues on stochastic models for the velocity of the fluid seen 


At the end of the analysis of stochastic models for polydisperse two-phase flows, it is worth 
emphasizing that modeling the velocity of the fluid seen remains an open issue. Coming 
up with tractable models that still represent proper descriptions of the physics of turbulent 
flows is not an easy task. Indeed, it must be remembered that the velocity of the fluid seen 
is Us(t) = U(f,Xp(f)) which means that, though Us(t) is a particle-attached variable, it 
nevertheless involves time and space correlations of the carrier turbulent fluid flow. Thus, 
the crossing-trajectory effect (GTE) { 4 ], which is related to particle and fluid velocity slips, 
induces several challenges for the formulation of a model for Us(f). 

With respect to this situation, the standpoint chosen in the present work is to consider 
jointly a set of criteria that points to acceptable forms of the different terms entering a 
stochastic diffusion model retained for the velocity of the fluid seen. As the complexity of 
the physics involved can, unfortunately, lead to overlooking the basic properties of a model, 
the criteria (P-1) and (P-2) have been useful to clarify some situations. Then, based on 
proper formulations of the fluid limit and on the GLM (which is the essence of (P-3)), (P-4) 
is important for the closure of the diffusion coefficient, (P-5) helps to bring out a correct 
form of the return-to-equilibrium term of a Langevin model, while (P-6) indicates possible 
expressions for the additional term of the drift vector related to the slip velocity between 
the fluid and the particles. Yet, they remain suggestions and indications. It must not be 
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forgotten that present models still contain assumptions. More precisely, the timescales of the 
velocity of the fluid seen are inputs (relying on the Csanady’s expressions) and, to the best 
of the authors’ knowledge, these expressions have not been worked out from hrst principles 
(see detailed discussions in sections 7.4 and 9.3.4 in Minier and Peirano jd] or hrst attempts 
in ^). This has consequences. For example, (P-6) is helpful to reveal the existence of a 
supplementary drift term, which basically comes from a hrst-order development of the huid 
velocity seen in the limit of small St. However, an important point is that it may not be 
relevant to push the analysis too far concentrating on that sole issue, when other aspects 
are disregarded. For example, arguing on whether the mean or instantaneous slip velocities 
should enter the additional drift term becomes secondary if the diffusion coefficient is not 
correctly closed and while it remains unclear whether the difference is accounted for (or not) 
by the Csanady’s formulas. The main messages are therefore: hrst, it is essential to address 
together the different terms entering a stochastic model and, second, that new ideas based 
on more fundamental derivations would be of value. 


VI. CONCLUSION 


In this paper, a new approach which consists in selecting a set of criteria has been de¬ 
scribed for Lagrangian stochastic models. Gathering present knowledge into a comprehensive 
set of guidelines has been shown to be useful both to assess whether existing models sat¬ 
isfy basic properties and also to help future developments. For that purpose, criteria have 
been put forward for Lagrangian stochastic models in single-phase (see (F-1) to (F-3) in 
section [Til Ap and in two-phase hows (see (P-1) to (P-6) in section ITlIBD . In the single-phase 
how situation, the emphasis is basically put on requiring that stochastic models be fully con¬ 
sistent with the high Reynolds-number structure of the Reynolds-stress equations. In the 
two-phase how situation, as the available information is diherent, the criteria insist on re¬ 
specting the particle-tracer limit, as well as basic properties of turbulence and what models 
should stand for. Although these criteria are somewhat diherent in the single- and two- 
phase cases, the approach followed in the present paper represents an attempt at addressing 
stochastic models devised for each situation from a unihed standpoint. 

It is believed that the present list of criteria is made up by simple and physically- 
meaningful requirements that must be met by Lagrangian stochastic models. Yet, from 
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the analyses which are summarized in Table [I] for single-phase models and in Table [TTl for 
two-phase models, hrst conclusions can be drawn. For single-phase flows, it was brought out 
that present normalized Langevin models are flawed since they do not respect convective 
transport in general non-homogeneous flows, as shown in section HV A 21 For two-phase flows, 
the situation is more confused and, in that respect, the analysis carried out in section IV Al is 
helpful to clarify the modeling picture. It appears that, among those considered there, only 
one formulation is acceptable, in the sense that all the criteria chosen in section IIII Bl are 
met. Not surprisingly, this model remains the only two-phase PDF model which has been 
validated in engineering conhgurations jb]. This is clearly a very poor and unsatisfactory 
state and much work remains to be done. In that sense, another conclusion is that a safe 
approach is to formulate models in terms of instantaneous fluid particle velocities, as for the 
GLM in single-phase flows. 

It is hoped that the guidelines put forward here will lead to improved Lagrangian stochas¬ 
tic models. For example, it appears that the complete Langevin model in two-phase flows 
is based on the SLM and, therefore, one option could be to devise extensions to obtain the 
counterpart for the velocity of the fluid seen of the GLM for the velocity of fluid particles. 
Furthermore, the present set of criteria is not meant as a dehnitive choice and could be 
extended in further works. For instance, low-Reynolds, Brownian effects or specihc aspects 
of particle preferential concentration effects could be considered through the formulation of 
new criteria to guide model assessment in these limits. 

Another important aim of the present approach is to provide guidelines for future model 
developments and in related subjects. In particular, the present analysis has been carried 
out in the framework of RANS approaches where a complete formalism is available. Yet, the 
methodology followed here, as well as the conclusions which have been reached, have direct 
implications for particle-laden turbulent flows simulated with a LES (Large-Eddy Simula¬ 
tion) method. Indeed, in the LES approach, a specihc model is needed to account for the 


unresolved part of the huid velocity seen by discrete particles 
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74| which corresponds 


to the huctuating velocity in the RANS context. Among the few models developed so far, 
most are written in terms of this unresolved part of the huid velocity seen [^, 75| but some 
formulations B include terms similar to the ones appeariiig on the rhs of Eq. fl39bl) while 
others retain only the simplest form of a Langevin equation 7^. Furthermore, some propos¬ 
als disregard ehects due to particle inertia or the GTE ehect and write an equation with a 
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le unresolved part 


diffusion coefficient as for fluid particles 0 whereas others assume that t 
is at equilibrium to close the diffusion coefficient of the Langevin model 7^. On the other 
hand, it is interesting to note that, for engineering applications, another proposal 0,0] was 
formulated in terms of the instantaneous fluid velocity seen with a model that is an extension 
to the LES context of the complete Langevin model discussed in section IV A 21 with corre¬ 
sponding expressions of the Csanady’s time scales and diffusion terms but without the mean 
slip term in the drift vector that is essential to the criterion (P-6). Clearly, there is some 
uncertainty and the issues addressed in this paper concerning the formulation of Lagrangian 
stochastic models (instantaneous versus fluctuating or unresolved fluid velocity components, 
closures of the drift and diffusion terms in the two-phase flow case, etc.) are relevant for the 
question on how to express subgrid effects to simulate two-phase flows with the LES method. 
However, as the statistical operator involved is different (a spatial hltering is applied instead 
of a mathematically well-defined probabilistic expectation), a rigorous formalism is needed 


to address these issues. The situation is more advanced for single-p 


rase flows where the FDF 


(Filtered Density Function) formalism has been been developed 
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79|. First steps have 


been proposed in the two-phase case 80|| but additional work is needed to clarify the for¬ 
mulations of Lagrangian stochastic models in the LES approach to particle-laden turbulent 
flows. 
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(F-1) 

(F-2) 

(F-3) 

Generalised Langevin model ^ 

/ 

/ 

/ 

Incomplete Langevin model for fluctuating velocities ^ 

/ 

X 

/ 

Incomplete Langevin model for fluctuating velocities ^ 

/ 

/ 

X 

Complete Langevin model for fluctuating velocities ^ 

/ 

/ 

/ 

Normalised Langevin model ^ 

X 

/ 

X 


^ as analysed in section HV Al 

as analysed in section HV A II where the term labelled (a) in Eq. (|39bl) is missing 

as analysed in section HV A II where the term labelled (b) in Eq. (|39bl) is missing 

as analysed in section IIV A II where all the terms on the rhs of Eq. (I39bl) are present 

® as analysed in section HV A 21 

TABLE I. Summary of the properties of the stochastic models for single-phase flows considered in 
section IIVI with respect to the criteria listed in section IIII A[ 



(P-1) 

(P-2) 

(P-3) 

(P-4) 

(P- 5) 

(P-6) 

(P-6bis) 

First Langevin model ^ 

/ 

X 

- 

- 

/ 

X 

/ 

Complete Langevin model (2001) ^ 

/ 

/ 

/ 

/ 

X 

/ 

/ 

Complete Langevin model (2004) ^ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

Normalised Langevin model ^ 

X 

/ 

X 

/ 

/ 

X 

X 

Hybrid-DNS model ^ 

/ 

X 

X 

X 

/ 

X 

/ 


as analysed in section FV A II (where the sign (-) indicates that the criterion cannot be checked) 
as analysed in section FV A 21 
as analysed in section FV A 21 
as analysed in section FV A 41 
® as analysed in section FV A 51 

TABLE 11. Summary of the properties of the stochastic models for two-phase flows considered in 
section |V] with respect to the criteria listed in section IIIIBI 
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